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ABSTRACT 

We have measured the column density distribution function, f(iVHi), at z = using 
21-cm Hi emission from galaxies selected from a blind Hi survey. f(iVm) is found to be 
smaller and flatter at z = than indicated by high-redshift measurements of Damped 
Lyman-a (DLA) systems, consistent with the predictions of hierarchical galaxy for- 
mation. The derived DLA number density per unit redshift, dN^LA/dz =0.058, is 
in moderate agreement with values calculated from low-redshift QSO absorption line 
studies. We use two different methods to determine the types of galaxies which con- 
tribute most to the DLA cross-section: comparing the power law slope of f(Ani) to 
theoretical predictions and analysing contributions to cLNdlai 'dz. We find that com- 
parison of the power law slope cannot rule out spiral discs as the dominant galaxy 
type responsible for DLA systems. Analysis of dN dla/ 'dz however, is much more dis- 
criminating. We find that galaxies with log A/hi < 9.0 make up 34% of dNoi,A/dz; 
Irregular and Magellanic types contribute 25%; galaxies with surface brightness, /I 25 > 
24 mag arcsec -2 account for 22% and sub-L» galaxies contribute 45% to dNoLA/dz. 
We conclude that a large range of galaxy types give rise to DLA systems, not just 
large spiral galaxies as previously speculated. 
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1 INTRODUCTION 

Lyman-a (Lya) absorption systems, in the continuum spec- 
tra of background QSOs, provide a powerful probe of the dis- 
tribution of neutral hydrogen in the Universe. Damped Lya 
(DLA) systems have column densities A?hi ^ 2 x 10 20 cm -2 
and contain the bulk of the ne utral gas in the Universe 
jStorrie-Lombardi fc Wolfell2000lh The value of Q. gas from 
DLAs at high redshift agrees well with Q s tar in present day 
galaxies. The fact that Q gas decreases monotonically with 
decreasing redshift and approaches Q ga s(z = 0) from 21- 
cm emission of local galaxies also provides strong ev idence 
that DLAs are caused by galaxies dWolfe et alj fl995). The 
properties of these galaxies however, remain mostly uncon- 
strained. 

There are two primary competing theories on the evo- 
lution of neutral gas in galaxi es. The first is mo nolithic col- 
lapse of protogalacti c discs JWolfe et alJll986l) . This sce- 
nario, supported by iProchaska fc Wolfd ((1997) , proposes 
that the kinematics of neutral gas tracing-ions in DLAs 
are best described by thick, rotating galactic discs. Alter- 
natively, neutral gas in galaxie s could be accre t ed hie rarchi- 
cally from sub-galactic halos. lHaehnelt et alJ (ll99STl argue 



that the kinematics of DLAs can equally be explained by 
sight-lines through irregular protogalactic clumps. 

Since massive spiral galaxies contain the bulk of the 
Hi at z = 0, it has been presumed that they also account 
for most of the A^hi ^ 2 x 10 20 cm ~ 2 cross-section (e.g. 
IWolfe et alJll986l : iRao fc Briggslll993Th With the advent of 
space-based UV spectroscopy, and the discovery of low red- 
shift [z < 1.65) DLAs, this assumption can now be tested. 
At low redshift the galaxies responsible for DLAs can be 
resolved, so their luminosities, surface brightness and mor- 
phological types can be determined. 

At z < 1.65 fewer than 30 DLAs are known. Images 
of 15 DLA galaxies have been published. The morphological 
breakdown includes 7 low surface brightness (LSB) galaxies, 
6 spirals and 2 compact dwarf-like galaxies <ILe Brun et al, 
| l997t|Pettini et alJ|2000 | ;|Turnshek et al.ll200lHBowen et~ 
2001; iNestor et alJ 2002). These galaxies are distributed 
evenly in luminosity, with 7 galaxies having L < L». Only 
five of these galaxies have measured spectroscopic redshifts. 
Unless the galaxy has a confirmed redshift matching that of 
the absorption systems, it is not certain that a particular 
galaxy is responsible for a DLA system. 

To determine the Hi emission properties of these 
DLA galaxies, they need to have even lower redshifts 
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(due to the weakness of the 21-cm transition). There 
is only one case of Hi emission detected from a DLA 
galax y. The DLA galax y SBS 1543+593 ijBowen et alJ 
| 200l|) has been detected l|Bowen et alJboOll) and mapped 
JChenealur fc Kanekarl2002^ in 21-cm. The candidate (Lya 
column dens ity yet to be p ublish ed) DLA galaxy NGC 4203 
reporte d bv IMil ler et al l lll999l) has also been mapped in 
21-cm Jvan Driel et al.lll988l) . 

For a set of random QSO sight-lines, DLA absorbers 
simply map out the cross-section on the sky with Am ^ 
2 x 10 20 cm -2 . For this reason, it is not surprising to find 
a variety of galaxies responsible for DLAs. There are also 
theoretical suggestions that different galaxy types, such as 
LSB galaxies may contribute significantly to the high Am 
cross-section Oimenez et allll999l : lBo"issier et al.l2 002l. One 
way to characterize this cross-section is to use the column 
density distribution function, /(Am). 

/(Am) is defined as the number of absorbers in the 
range Am to A^m + dNm per unit Am per unit absorption 
distance. The most recent determination of /(Am) over a 
range of redshifts, using a statis tical sample of 6 46 QSOs, 
with 85 measured DLAs, is by IStorrie-Lombardi fc Wolfd 
(2000). At low redshift, DLAs are rare due to the low vol- 
ume of space surveyed by QSO sight-lines and the fact the 
spectral observations must be made in t he UV . To increase 
the number statistics, IRao fc Turnshekl J200CT) used a alter- 
native approach by selecting candidate DLAs from a sample 
of Mgll absorption line systems. Although < 30 DLAs are 
now known at z < 1.65, the number at z < 0.1 is still only 
two. Significant evolution of galaxies (and hence DLAs) is 
still expected in the range < z < 1.65. Therefore DLAs 
are unlikely to provide a suitable technique for measuring 
/(Am) at z = 0. 

Fortunately, at z — the alternative method of 21- 
cm emission can be used to measure /(Am). The 21-cm 
emission method has two main advantages. Unlike QSO ab- 
sorption observations, 21-cm emission does not suffer from 
the low number statistics. Secondly, since the properties of 
galaxies detected in 21-cm are known, they can be used to 
infer the distribution of DLA galaxy properties such as lu- 
minosity, surface brightness, and Hi mass. Hence we are able 
to assign the probability that a particular galaxy type gives 
rise to a DLA at z = 0. 

Previous calculations of the column density distribu- 
tion function f rom 2 1-cm emission have been published by 
IRao fc Briggsl lll993Fl and IZwaan et alJ \ 19991. l2002Tl . The 
study bv lRao fc Briggsl fe^j) used Arecibo obse rvations of 
27 large spiral galaxies. IZwaan et al] (^999, 2002}) used syn- 
thesis observations of optically-selected Ursa Major cluster 
members and improved significantly on the previous calcu- 
lation. The benefit of using a galaxy cluster is that the spa- 
tial resolution is the same for each galaxy. H owever, galaxy 
clust ers are known to be Hi-deficient (e.g. Sola nes et al] 
2001) and may also have fl atter Hi mass functions (e.g. 
iRosenberg fc Schneiderl l2002l which could bias the /(Am) 
calculation. 

To calculate /(Am) from 21-cm emission we must 
be confident that 21-cm measurements of Am are 
equivalent to those from Lya absorption. According to 
iDickev fc Lockmanl dl99ff) "The agreement between ultra- 
violet and radio estimates of Nm is excellent and somewhat 
astonishing, considering that the angular resolution, experi- 



mental techniques, and telescopes differ in all respects" . This 
statement refers to a comparison of high latitude Galac- 
tic 21-cm measurements compared with Lya absorption to- 
ward distant stars (in the same direction) in the range 
10 20 < Nm < 10 21 cm~ 2 . A further issue is that high col- 
umn density gas may clump on scales less t han the beam 
size. This concern has been addressed by IRao fc Briggsl 
lll993l) in their calculation of the contribution of high Am 
gas to the Hi mass function (which is not affected by beam 
size) compared with the integral of mass from the /(Am) 
cross-section. They find the integral of mass in /(Am) is 
equivalent to that from the mass function at high AThi ■ They 
concluded that any high Am column density gas missed due 
to a large beam size is negligible. 

In this paper we calculate /(Am) from a randomly se- 
lected Hi sample of galaxies selected from the whole south- 
ern sky. This selection process most closely mimics that of 
random QSO sight-lines. The selection process and observa- 
tions are described in §2. In §3 /(Am) at z = is calculated 
and compared to higher redshift measurements. The data 
set is in then sorted by galaxy properties and the contri- 
bution to the cross-section of Hi gas satisfying the DLA 
criterion is calculated. In §4 the implications of the results 
are discussed and conclusions are presented in §5. Details of 
the galaxy sample, including Hi contour maps and column 
density histograms are given in the appendix. 



2 OBSERVATIONS 

Galaxies for this study were selected from the Hi Parkes 
All-Sky Survey (HIPASS) Bright Galaxy Catalogue (BGC). 
HIPASS is a blind survey for Hi covering the entire south- 
ern sk y. The data reduction is described by Barnes et al. 
J2001]) . The 3cr Hi mass limit of HIPASS is ~ 10 6 D 2 Mpc M e , 
the 3er Hi column density limit is 4xl0 18 cm -2 and the 
FWHP beamwidth is 15.5'. The BGC is the largest blind Hi 
cxtragalactic catalogue to-date. The galaxy find ing ja ,nd pa- 
rameterisation is described by Koribalski et al. (1200211 . The 
Hi masses in the BGC were calculated using Local Group 
corrected velocities and an Ho = 75 km s _1 Mpc -1 . 

HIPASS poorly resolves most galaxies and can usually 
only measure a lower limit to the column density, equivalent 
to the Hi being smoothly spread over the entire beam. To 
measure the distribution of Hi column densities, higher reso- 
lution mapping is required. Thirty-five galaxies were chosen 
at random from the BGC for mapping with the Australia 
Telescope Compact Array (ATCA) 1 . To provide a represen- 
tative range of Hi mass, galaxies were chosen with an even 
mass distribution in the range 7.5 < log Mm< 10.5 M©. 
Galaxies were restricted to those with heliocentric velocities 
in the range 500 to 1700 km s" 1 so that the spatial resolu- 
tion remained comparable across the sample. 

Each galaxy was observed with the ATCA 375 and 750D 
array configurations. The data set was reduced in MIRIAD, 
the velocity resolution in each datacube is 3.3 km s _1 . The 
datacubes have an average restored beam of l'xl'and typ- 
ical RMS noise in line-free regions of 3.7 mjy beam -1 , this 

1 The Australia Telescope Compact Array is part of the Australia 
Telescope which is funded by the Commonwealth of Australia for 
operation as a National Facility managed by CSIRO. 
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corresponds to a 3er column density limit, over 3 velocity 
channels, of 2xl0 19 cm -2 . The Hi column density (iVm) of 
each pixel in each galaxy was calculated by summing flux 
density along the spectral axis. The integrated flux density, 
Sdv (Jy km s" 1 ), was converted to a brightness tempera- 
ture in the usual manner, T b = 1.36 x 10 3 A 2 S/b ma] b mm , 
where A is in cm and the major and minor beam axes, 
bmaj and b m in, in arcsec. Tb is converted to column den- 
sity using N m = 1.823 x 10 18 / T b dv cm " 2 . This calculation 
assum es that the gas is optically thin. iDickev fc Lockmanl 
( 1990) used typical Galactic Hi parameters (T s ~ 50K 
and Av ~ 10km s" 1 ) to show that the optical depth is 
greater than 1 for Am > 10 21 cm -2 . Consequently the 
column density may be underestimated as Am approaches 
10 21 cm -2 , providing an upper limit to these observations. 
Regions of high Am may also be underestimated due to 
the resolution of the observations. Nearby galaxies which 
have been mapped at high resolutio n show many regions 
with A m > 10 21 cm" 2 (e.g. the LMC. IStavelev-Smith et alJ 
2003). Table Al summarises the Hi parameters for each of 
these galaxies. Also included are optical parameters from 
the Lyon-Meudon Extragalactic Database (LED A). Some 
galaxies are new or recent discoveries and therefore some 
optical parameters are unavailable. For example, HIPASS 
J0949-56 lies behind the Milky Way and no optical counter- 
part is visible due to high foreground extinction (As=9.2). 
Figure lAt gives the Nm contours overlaid on a DSS image 
of each galaxy. An Am histogram for each galaxy is given 
in Figure lA"fe. 



3 DATA ANALYSIS 

The column density distribution function, from each galaxy 
(i) in the sample, is calculated using 
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where A(7Vhi) is the area in Mpc 2 occupied by column den- 
sities in the range -Am to Am + dNm (d logAm = 0.2 dex is 
used throughout). The function is normalised by the space 
density, <j>(Mm) t aken f rom the HIPASS BGC Hi mass func- 
tion (Zwaan et al. 120031) . which has a faint end slope of a = 
— 1.30±0.08 and a characteristic Hi mass of log(A'/*/Mo) = 
9.79±0.06 and normalisation 6* = (8.6±2.1) x 10~ 3 Mpc" 3 . 
The resulting distribution function is plotted in Figure A 
power law, /(Am) = BAm"' 3 , is fitted with (3 = 1.5 ± 0.1 
in the range 19.6 < logAm < 21.6. The three highest 
column density points do not agree well with this power 
law. A double power law is therefore fitted with a break 
at log Am = 20.9. In this case the lower column densities, 
19.6 < logAm < 20.9, are best fitted with (3 = 1.4 ± 0.2 
and the upper column densities, 20.9 < logAm < 21.6, 
with (3 = 2.1 ± 0.9. The double power law is also given in 
Figure Q 

The z = /(Am) from 21-cm is compared to higher 
redshift measurements from QSO absorption line studies. 
The high r edshift ((z) ~ 2.5) DLA meas urements are 
taken from IStorrie-Lombardi fc Woifel <|2000T) and the in- 
termediate redsh i ft ((z) ~ 0.78) DLA measurements from 
iRao fc Turnshekl j200Cf) . (Although qo is unimportant at 




, <z> = (This work) 
, <z> = 0.73 (Rao & Turnshek, 2000) 
,<z> = 2.5 (Storrie-Lombardi & Wolfe, 2000) 



bo, HI column density (cm 2 ) 

Figure 1. Column density distribution function at z = for all 
galaxies in the sample (data points). The error bars are calculated 
using Poisson statistics. A comp arison is made to the high redshift 
({z) = 2.5) measurements of IStorrie-Lombardi fc Woifel (2000) 
(t riangles) and intermedi ate redshift ((z) = 0.78) measurements 
o flRao fc Turnshekl lioOO'l (circles). The distribution is fitted with 
a single power law (solid line) and double power law (dotted line). 



z = 0, the high redshift /(Am) and dNDLA/dz use qo = 0.) 
We find that with decreasing redshift, /(Am) decreases at 
all column densities and flattens for 20.9 < log-Am < 21.6 
from (z) ~ 2.5 to 0. Quantitatively, the flattening is deter- 
mined by fitting a power law in this high Am range to the 
(z) ~ 2.5 (P = 2.5) and z = (f3 = 2.1) distribution func- 
tions. The flattening of /(Am) at high column densities is 
attributed to growth of denser objects (i.e. galaxies), there- 
fore reducing the number of lower column density systems 
and increasing the high Am systems. The overall decrease in 
/(Am) to z = can be explained by the c onsumption of Hi 
gas to form stars. The magnitude of the IRao fc Turnshekl 
(2000) /(Am) fits with the general evolutionary picture. 
However the slope of the function is not in good agreement, 
although the error bars are large enough for the above de- 
scription to hold within the stated uncertainties. 

The departure from the power law at .Am < 10 20 cm~ 2 
is perhaps partly due to the H i-Hh transition , descr ibed as 
the "footprint" of ionisation bv lCorbelli et all <l200ll) . essen- 
tially marking the Hi edge of the galaxies. It may also be 
due to the detection limit of our observations. The detection 
limit is a function of the line width, so measurements of low 
column density regions of galaxies with large line widths are 
most likely underestimated. 

Systematic uncertainties in /(Am) could arise from er- 
rors in the galaxy distances or the Hi mass function. The 
distance to each galaxy was calculated using the radial ve- 
locity, corrected for the motion of the Local Group. Alter- 
natively, a parametric model for the velocity field of galax- 
ies cajn_l2e_jised tocalculate distance. The model used is 
bv lTonrv et al] <l200Cf) . uncertainties in the distance affects 
the </>(M H i) and A(Nni) terms in /(Am). The difference 
between the two distance calculations is ~ 20% for most 
galaxies in the sample. The resulting /(-Am) is plotted in 
Figure along with the original calculation. In addition, 
/(Am) is calculated wit h two different Hi Mass functions, 
one with shallower slope dZwaan et alJll997Tl and one with a 
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Figure 2. The effect of systematic errors on the column den- Figure 3. Column density distribution function for individual 

sity distribution function. The original fu nction is included spiral galaxies only (solid lines). The dotted line, Nm~ 3 , is pre- 

along with two different Hi mass functions IZwaan et al Jtl997t dieted theoretically (arbitrary normalisation). 

IRosenberg fc SchnpiHfiti2002^ a nd a velocity flow m odel to calcu- 
late the distance to each galaxy iTonrv et al.l200cl) . The different 
functions show that systematic uncertainties are small and have 
little effect on the shape of /(-/Vjji)- 



stepper slope (IRosenberg fc Schneiderl2 002). These system- 
atic uncertainties have little effect on the shape and normal- 
isation of the column density distribution function beyond 
the random errors already quoted. 

The power law slope of the column density distribu- 
tion is theoretically predicted to be —3 at the hig h Nm 
end for disc galaxies (lMilgromlll988l : IWorfe et alJl 995). This 
prediction together with a measured f(Nm) slope natter 
than —3 has been used to argue t hat disc galaxies canno t 
be responsible for all DLA systems (iRao fc Turnshekfe oOO ) . 
Considering only the spiral galaxies in the sample, we find 
that f(iVm), for Nm^? 10 21 cm -2 , has a power law slope of 
—2.6 ± TO. In Figure [3J f(-/Vm) for each spiral galaxy con- 
tributing to this Nm range is shown. If we consider both 
spiral and non-spiral galaxy types, the power law slope in 
this high iVm region increases to —2.1 ± 0.9. Thus, includ- 
ing non-spiral galaxies in the sample, the data still agrees, 
within the uncertainties, with the theoretical prediction for 
spiral galaxies only. 

To emphasise the relationship between f(Nm) and to- 
tal galaxy Hi mass these are plotted in a range of column 
density bins in Figure 2] f(Nm, Mm) has the same basic 
shape at all column densities. For most Nm, f(Nm, Mm) is 
essentially a flat function of Hi mass, rising slightly from log 
Mm = 9 to 9.5 Mq and dropping off sharply after log Mm 
= 10 Mq. To compare with DLA absorption line statistics, 
we also integrate f(Nm) from TVhi ) 2 x 10 20 cm -2 , which 
is equivalent to calculating the cross-section of Hi satisfying 
the DLA criteria, or the number density of DLAs per unit 
redshift. 

^pA = f Nm - f{Nm)dNm (2) 
dz JiV HT= 2xio 2 o 

For the whole sample, cLNdla/ 'dz = 0.058T0.006 (random) 
-o 02 (systematic) at z = 0. The random uncertainty is 
based on Poisson statistics and the systematic uncertainty 



is from the Hi mass function. This value agrees very well 
with the parameterisation , dNpLAl 'dz= 0.055(l-|-z) 1 1 1 over 
a large redshift range bv IStorrie-Lombardi fc Wolfd ll200fj) 
(see their Figure 11). This r esults confirms the suggestion by 
IStorrie-Lombardi fc Wolfel (|200(t l that there is no intrinsic 
evolution in the product of space density and cross-section 
of DLAs with redshift to z = 0. 

Figure compares the fractional contribution of 
dNoLA/dz to the Hi mass density, pm=0(Mm) x Mhi as 
a function of log Mm- The histogram shows that 34% of 
dNoLA/dz can be attributed to galaxies with log Mm < 9-0. 
This is in contrast to pm (continuous line in figure), which 
is dominated by M* galaxies. In comparison 26% of pm is 
due to galaxies with log Mm< 9.0. Furthermore, Figure H2 
shows that the mean DLA cross-section of log Mm < 9.0 
galaxies is a mere 4%. 

Although the cross-section of Hi satisfying the DLA cri- 
terion in an individual low Hi mass galaxy may be small, 
the product of cross section and space density is more sig- 
nificant than the product of Hi mass and space density. This 
leads to dNoLA/dz as a function of Hi mass peaking at less 
than M*. Low mass galaxies may contribute little to the Hi 
mass density, however they make up an important part of 
dNoLA/dz. 

This result is similar to Ro senberg fc Schneider! (120031) 
if the different Hi mass functions are considered. In both 
cases the dNoLA/dz histogram generally sits slightly above 
Phi for Mm < M* and slig htly below for Mm > M* . 
IRosenberg fc Schneider! d2003l) find a tight correlation be- 
tween log Mm and log A(Nm > 2 x 10 20 cm -2 ) and use 
this to calculate dNoLA/dz. In that case pm oc dNoLA/ 'dz, 
so the shape of pm(Mm) should match dNoLA/dz{Mm) 
peaking at M* . In practice the shape of dNpLA/d z devi - 
ates slightly from pm since IRosenberg fc Schneider! j2003ll 
use the 1/Vtat value for each individual galaxy rather than 
4>(Mm) from the function form of the Hi mass function. 



The Column Density Distribution Function at z 



= from Hi Selected Galaxies. 



5 



• 9.5 
9 9 



0.8 



neon DLA cross section 



7.0 7.5 8.0 8.5 9.0 9.5 10.0 10.5 

log (m h / m ) 

Figure 4. Column density distribution function at z = as a 
function of galaxy Hi mass, plotted at six different column den- 
sities. 
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Figure 6. Fractional contribution of each Hi mass bin to 
dNpLA/dz (solid line) and the mean DLA cross-section (dashed 
line). 
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Figure 5. Fractional contribution of each Hi mass bin to 
dNuLAl 'dz (histogram) and Hi mass density, pm (Mir,Mpc~ 3 , 
continuous line) comparing the work in this paper (solid his- 
togram and line) with Rosenberg & Schneider (2003) (dashed 
histogram and line). 



3.1 Relationship between DLA-gas and Optical 
Properties 

In this section we explore the relationship between 
dNoLA/dz and optical properties of galaxies in the sam- 
ple. dNoLAl 'dz is still normalised by the Hi mass function, 
but is weighted to reflect the number distribution of opti- 
cal properties in the sample. The mean DLA cross-section, 
A(Nm > 2 x 10 20 cm" 2 ), for galaxies with each different 
property is also calculated. The fractional contributions to 
dNoLA/dz and mean DLA cross-section as a function of 
each galaxy property is given in Figures Q to I1UI In each 
case the dNoLA/dz contribution (solid histogram) describes 
the probability of a galaxy with a particular property giving 
rise to a DLA system. The mean DLA cross-section (dashed 
histogram) indicates for a galaxy with a particular property, 
its non-weighted DLA cross-section on the sky. 

First, the likelihood of a galaxy in different luminosity 



ranges giving rise to a DLA system is calculated. Galaxies 
in the range 8 Lq < log Lb < 10 Lq dominate dNoLA/ 'dz. 
Sub-L» (log Lb < 9 Lq) galaxies account for 45% of the 
DLA cross-section. In contrast these sub-L* galaxies have 
very small (6%) DLA cross-sections compared with L > L, 
galaxies. This has interesting consequences for the interpre- 
tation of DLA-galaxy associations. For example, although 
sub-L, galaxies are just as likely to give rise to a DLA ab- 
sorber, given a DLA system in a field with one sub-L, and 
one super-L* galaxy, the super-L, galaxy is ~10 times more 
likely to be responsible for the DLA by virtue of its much 
larger Hi cross-section. 

To investigate the relationship between dNoLA/dz and 
galaxy morphology, the sample was divided into three 
parts. The first group contains Irregular and Magellanic- 
type galaxies. The second group, late type spirals, includes 
Sd and Sc types. The early type spiral group is comprised 
of Sb, Sa and SO types. We find that dNoLA/ 'dz is dom- 
inated by late type spirals, which contribute 48% to the 
cross-section of DLA-gas on the sky. The contribution of 
Irregular and Magellanic types is non-negligible at 25%. In- 
dividually however, early type spirals each have the largest 
DLA cross-section (40%). 

dNoLA/ 'dz as a function of galaxy mean surface bright- 
ness within the magnitude 25 isophote (7*25 ) i s dominated 
by galaxies with 23 < Jl 25 < 24 mag arcsec -2 . The mean 
DLA cross-section of these galaxies also dominate the sur- 
face brightness division. Lower surface brightness galaxies, 
with Jl 2 5> 24 mag arcsec -2 , make up 22% of the cross- 
section. Five galaxies did not have this parameter available 
in the LEDA database. One of these is HIPASS J0949-56, 
which lies behind the Milky Way. Another, HIPASS J1320- 
21 is a high surface brightness SO galaxy. The other three 
(HIPASS J0905-56, J1321-31 and J2222-48) all appear to be 
LSB galaxies (see Figure Al). The fact that these galaxies 
are missing from the analysis introduces a bias against low 
surface brightness galaxies. Including these three galaxies 
in the lowest surface brightness bin (25 < Jl 25 < 26 mag 
arcsec -2 ) increases the dNDLA/dz contribution from 3 to 
6 ± 5%. 

The relationship between dNoLA/dz, Hi mass and Lb 
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Figure 7. Fractional contribution of each luminosity decade to 
dNuLAI 'dz (solid line) and the mean DLA cross-section (dashed 
line). 
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Figure 9. Fractional contribution of each mean surface bright- 
ness decade to dN^LA/dz (solid line) and the mean DLA cross- 
section (dashed line). 
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Figure 8. Fractional contribution of each morphological type to 
dNoLA/dz (solid line) and the mean DLA cross-section (dashed 
line). 
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Figure 10. Fractional contribution of each mass-to-light ratio 
decade to dNjjLA/dz (solid line) and the mean DLA cross-section 
(dashed line). 



is combined in the calculation of Mm-to-La ratio. The cross 
section of DLA-gas is dominated by gas-rich galaxies with 
Mm-to-Ls ratios in the range 1 - 10. The very high Mm-to- 
Ls ratio galaxies contribute little to the DLA cross section 
(5%), mostly due to their rarity. It is also due to the fact that 
most of the Hi in these galaxies has a low column density. 
The contribution of very high Mm-to-La ratio galaxies to 
dN/dz for log A>m sC 19.1 cm" 2 is 18 ± 5%. 



4 DISCUSSION 

The flattening of f(A*Hi) between (z) — 2.5 and z = for 
Am > 10 20 ' 5 cm -2 can be explained by accreting clumps 
of Hi into the central regions of galaxies and DLA sys- 
tems. Milgrom (1988) argues that the slope of the column 
density distribution function reflects the density distribu- 
tion within the absorbers themselves. A case is presented 
where lines of sight through spherical gas clouds give rise to 
i(Nm)oc Am - ' 3 . The column density profile of each sphere 



is given by Nm(l) = iV [l + (l/r c ) 2 ]~ c , where I is the im- 
pact parameter and r c is the core radius. The power law in- 
dices of the local and global distribution are related through 
/3 = (c+ l)/c. The high redshift c olumn density distribution 
JStorrie-Lombardi fe Wolfel l200d) . for A^hi ^ 10 20 ' 5 cmT 2 , 
has a power law index of /3 = 2.2. Comparing /3 = 2.2 (with 
c = 0.8) to our z = measurement of (5 — 1.6 (with c = 1.7), 
yields a more centrally concentrated column density profile 
for galaxies (or DLAs) at z = 0. This picture lends support 
to the hierarchical construction of galaxies and DLAs. The 
alternative model, monolithic collapse into disc-like systems, 
predicts f(Nm) oc A^hi -3 for very high column densities 
(Ahi^ 10 21 cm" 2 ). Our data is marginally consistent with 
this theoretical prediction for disc galaxies. However, we are 
unable to confirm a shallower power law for non-disc galax- 
ies, because of a lack of measured high column densities, due 
to resolution limitations. 

The number density of DLA systems per unit redshift 
is calculated to be dNoLA/dz — 0.058±0.006 (random) 
-o 02 (systematic) at z = 0. This result agrees well with 
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an extrapolation of higher redsh ift evolution measured by 
IStorrie-Lombardi fe Wolfel §000) . Our result is in moderate 
agreem ent (within the large uncertainties) with IChurchiTll 
fcOOll) who finds dN DLA /dz = 0.08±g;^ at (z) = 0.06 
from Mgll QSO absorption systems. Our result i s also 
con sistent with the 21-cm studies by IZwaan et alJ l|2002Tl 
and iRosenberg fc Schneiderl J2003f) who find dNoLA/dz = 
0.042±0.015 and 0.053±0.013 respectively. 

The association of Lya absorption lines with particular 
galaxies is a subject of much debate. Some statistical studies 
compare galaxy survey data complete to a limit of L* with 
Lya absorption systems (e.g. IPenton et alJ l2002:l . We have 
shown statistically, that sub-L* galaxies account for 45% of 
the DLA cross-section at z = 0. If such a large proportion 
of galaxies are missed, the strength of the galaxy- absorber 
cross-correlation function may be underestimated. 

The results presented here agree with the targeted 
searches for low redshift DLA galaxies that find a range of 
morphologies and luminosities. Unfortunately 21-cm emis- 
sion measurements of galaxies are limited to 2 < 0.1, 
whereas the detection Lya absorption is bounded only by 
the brightness of the background probe. However, we can 
project these results to redshifts beyond zero. We have 
shown that low Hi mass, low luminosity galaxies contribute 
significantly to the DLA cross-section at z = 0. These types 
of galaxies are more common at higher redshifts (e.g. lEllia 
1997). Provided the DLA cross-section of these types of 
galaxies do not get significantly smaller with redshift, it is 
likely they also account for an important part of the high 
redshift DLA cross-section. 

In future studies, where the absorption-line statistics 
and properties of z — DLA galaxies are better known, 
f(iVm) and dNoLA/dz from 21-cm emission will provide a 
test of possible biases in the DLA galaxy population. Se- 
lection effects include gravitational lensing of QSO-galaxy 
alignments, bringing otherwise dimmer QSOs into a statis- 
tical sample and creating a 'by-pa ss' effect where the line- 
of-sight avoids the galaxy cen tre iSmette et al.lll997Tl . Al- 
though, iLe Brun et alJ J200(J) showed that lensing is not 
important in the s ample of DLA galaxies presented in 
ILe Brun et alJ <ll997l) . they note that brighter QSOs may be 
affected. Spiral galaxies could be under-represented in the 
DLA galaxy population due to their large dust content, ob- 
scuring any background QSO JOstoiker & Heislcr 1984). A 
recent study of radio loud QSOs i Ellison et al.l l2001) found 
that a dust-induced bias in optical QSO surveys may have 
led to an underestimate of Qdla and dNoLA/dz by at most 
a factor of two. However, DLAs are 4 times a s likely to be 
found in the proximity of radio loud QSOs dEllison et alJ 
2002) which could overestimate Q.dla and dNoLA/dz. If 
the DLA cross-section from both 21-cm and QSO absorption 
is well measured and compared these biases can be better 
understood. 



5 CONCLUSIONS 

We have calculated the column density distribution function 
at z = from an Hl-selected sample of galaxies. We find 
that f(Nm) follows a power law with slope = —1.5 ± 0.1 
in the range 19.6 < logA^Hi < 21.6. When compared to 
high redshift QSO absorption line data, f(Nm) flattens and 



decreases with evolution to z — 0. This evolution can be in- 
terpreted as accretion of lower column density systems onto 
higher ones, thus reducing the number of low column den- 
sity galaxies. This picture fits with the hierarchical forma- 
tion model of galaxies and DLA systems. Conversely, we 
are in marginal agreement with the theoretical prediction 
f{N m ) oc N m ~ 3 for very high column densities in spiral 
discs. The fact that we find f(Nm) oc N m ~ 2 ' 1±0 ' 9 for very 
high column densities, for all types of galaxies, including 
spiral discs means that we cannot use this analysis to rule 
out the alternative model, that spiral discs are responsible 
for all DLA systems. 

However, analyzing the the slope of the f(Nm) power 
law is not the only method of determining the likely galaxy 
types which gives rise to DLA systems. The calculation of 
dNoLA/dz as function of Hi mass, luminosity, galaxy mor- 
phology, Mm-to-I/fl ratio and surface brightness reveals a 
range of galaxy types contribute to the cross-section of Ahi 
satisfying the DLA criteria. In particular we find that galax- 
ies with log Mm < 9.0 make up 34% of dNDLA/dz; Irregu- 
lar and Magellanic types contribute 25%; galaxies with sur- 
face brightness, Jl 25 > 24 mag arcsec -2 account for 22% and 
sub-L, galaxies contribute 45% to dNoLA/ 'dz. These results 
agree with findings from the imaging of low redshift DLA 
galaxies. Integrating over all galaxies in the sample we find 
that dNoLA/dz— 0.058, which is consistent with no intrinsic 
evolution in the product of space density and cross-section 
of DLAs with redshift. This result agrees with the redshift 
evolut ion predicted from high and int ermediate DLA sys- 
tems jStorrie-Lombardi fc Wolfel 12000) and is moderately 
con sistent with very low redshift QSO absorption line stud- 
ies (IChurchilll2o'o"lT> . 
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APPENDIX A: GALAXY SAMPLE 

Includes a table of all the galaxies, an optical image with 
Am contours overlaid and a column density histogram for 
each galaxy. 

This paper has been typeset from a TffX/ P/Tp^X file prepared 
by the author. 
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Table Al. Hi parameters for each galaxy from HIPASS and optical parameters from LEDA. The columns are the HIPASS 
name, optical ID, heliocentric velocity, distance (calculated using Local Group corrected velocities), Hi mass, beam size 
of 21-cm map with the minor axis in kpc in parentheses, mean surface brightness within the 25th magnitude isophotc, 
Afm-to-Lfl ratio and galaxy morphology. 



a Henning ct al. (?) 
b Banks ct al. (?) 
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Figure A2. Whi histograms for all galaxies in the sample. In each case, the x axis is log Nm (cm 2 ) and the y axis is the number of 

pixels. 



